 result=importdata("~/bse2020-3-8/bsegit/bse/search.out");
% csvwrite("~/bse2020-3-8/bsegit/bse/search.csv",result);

% 11 mx 超巨星质量 14 现在的轨道周期
% result1=result(result(:,3)>1000,:);
% result2=result(result(:,3)>100,:);
% result3=result(result(:,3)>1,:);
% result4=result(result(:,3)>1,:);
% plot(result4(:,12),log10(result4(:,14)),'b.')
% hold on
% plot(result3(:,12),log10(result3(:,14)),'g.')
% plot(result2(:,12),log10(result2(:,14)),'y.')
% a=plot(result1(:,12),log10(result1(:,14)),'r.');

% hold off
% apatch=patch(result(:,12),log10(result(:,14)),result(:,3),'edgecolor','flat','facecolor','none');
% apatch.LineStyle=   'none';
% apatch.Marker='.';
% scatter(result(:,12),log10(result(:,14)),3,log10(-result(:,end)),'filled')
xlswrite("~/bse2020-3-8/bsegit/bse/search20200325beta1.xls",result)
% scatter(result(:,1),log10(result(:,2)),10,log10(result(:,3)),'filled')
% xlabel('伴星质量/太阳质量')
% ylabel('log10(轨道周期)/log10(day)')
% colorbar
clear